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Abstract 

Using atomistic pseudopotential wave functions we calculate the electron 
and hole charging energies of InAs quantum dots. We find that the charging 
energies depend strongly on the dielectric constant eout of the surrounding 
material, and that when the latter is smaller than the dielectric constant of the 
dot (weak external screening) the electron-electron and hole-hole interactions 
are dominated by surface polarization effects. We predict the addition energies 
and the quasi-particle gap as a function of size and eout- We find excellent 
agreement with recent single-dot tunneling spectroscopy data for eout = 6. 
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Semiconductor quantum dots can be made with various dielectric coatings: Organic 
molecules other semiconductors (e.g. self-assembled dots [@, core-shell nanocrystals 

lithographically-etched dots 0, strain- induced dots 0), or glasses It has been 
realized that the dielectric environment can profoundly affect the optical and transport 
properties of quantum dots. This can be seen by considering the two processes described in 
Fig. 1, where a quantum dot of dielectric constant is embedded in a material of dielectric 
constant eout- Figure 1(a) depicts the process of adding three electrons to an otherwise 
neutral quantum dot. The initial configuration of the system, of energy Eq, consists of a 
neutral dot in the ground state and a Fermi reservoir at the reference energy e^ef = 0. The 
"charging energy" fii required to load the first electron into the quantum dot is 

= El — Eq = Eel + Sel\ (1) 

where Ei is the total energy of the dot with one additional electron, £ei is the energy of the 
single-particle level el with respect to the reference energy £ref, and Eg°^ is the self-energy 
of the additional electron interacting with its own image charge created by the dielectric 
mismatch at the surface of the dot 0. The charging energy /i2 to add the second electron 
to the quantum dot is 

fl2 = E2 — El = Eel + Sg°' + Jel,el , (2) 

where Jei,ei is the Coulomb interaction between the two electrons. It includes a direct 
electron-electron contribution J^i^l and a polarization contribution J^i^i arising from the 
interaction of one electron with the image charge of the other electron 0. Finally, the 
charging energy for the third electron is 

fJ'S = E3 — E2 = Ee2 + Se2^ + 2Jel,e2 — -^el,e2 , (3) 

where -^"61,62 is the exchange energy between two electrons with parallel spins in the el and 
e2 single-particle levels. The "addition energies" for the second and the third electrons are, 
respectively 

Ai,2 = /U2 - /il = Jel,el , (4) 
^2,3 = l^Z- 1^2 = (£e2 - £el) + iX'lt " ^if) + 

(2Jel,e2 ~ <^el,el) ~ -^el,e2 • (5) 

Since and J?°^ depend strongly on the dielectric constant of the surrounding material, 
the charging spectroscopy PJl^l of a quantum dot depends on its dielectric environment. 



Figure 1(b) describes the process of removing an electron from the highest occupied 
orbital of a neutral quantum dot and placing it into the lowest unoccupied orbital of an 
identical dot (located at infinite distance from the first dot). The energy required by this 
process ( "quasi-particle gap") is the difference between the ionization potential and the 
electron affinity of the dot. The initial configuration, consisting of the two neutral dots in 
the ground state, has energy 2Eq, while the final configuration has energy Ei + E_i, where 
E_i is the energy of the quantum dot with a hole in the highest occupied orbital hi. The 
quasi-particle gap is then 
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= E^ + E^,-2Eo = (Sei - e^) + S^f + S^f , (6) 

where £gap = ^ei — ^hi is the single-particle (HOMO-LUMO) gap. We see that the quasi- 
particle gap depends, via the polarization self-energies Sg°' and on the dielectric envi- 
ronment. The optical gap differs from the quasi-particle gap by the electron- hole interaction 

Jhl,el- 
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Very recently single-dot tunneling spectroscopy was applied to InAs nanocrystals , mea- 
suring An,n+i and e^Pp. 

The effects of dielectric confinement on the excitonic gap and the charging energies of 
quantum dots have been addressed in the past P,pJ]- |l5|] using the effective-mass approxi- 



mation (EMA). Recent pseudopotential calculations |TB| have demonstrated the importance 
of using an atomistic description of the quantum dot electronic structure for calculating the 
electron-hole Coulomb energy. The pseudopotential approach provides an accurate descrip- 
tion of the wave function decay outside the quantum dot and of the interband coupling due 
to quantum confinement, which are critical for a correct evaluation of the polarization and 
Coulomb energies in small nanocrystals. 

Using pseudopotential wave functions, we discuss here the effects of dielectric mismatch 
(cout 7^ Cin) on (a) the electron and hole charging energies fiN = En — E^-i, where E^ 
is the ground-state energy of electrons (or holes) in the quantum dot, (b) the addition 
energies A7v,Ar+i = fJ^N+i — I^-n, and (c) the quasi-particle band gap e^^p- We find excellent 
agreement with recent experimental results |10|], and interpret the data in terms of Coulomb 
and polarization contributions [Eqs. (|1|) - (H)]. We show that depending on the ratio between 
ein and Cout one encounters two physically distinct regimes of transport behavior: 

(i) When Cout -C ein (weak external screening) the electron-electron interaction Jtj is 
dominated by the polarization contribution JP°\ which depends only weakly on the single- 
particle states i and j. The charging energies /in depend strongly on the dielectric constant 
eout, and are widely spaced for various (large addition energies A7v,Ar+i). 

(ii) When eout > (strong external screening) the dominant contribution to Jij is the 
Coulomb energy J^j^^, which is quite sensitive to the identity of the states i and j. The 
charging energies depend weakly on eout and are more closely spaced (small addition 
energies) . 

The practical significance of these results stems from the fact that, due to the long- 
range character of the Coulomb interaction and the exponential decay of the wave functions 
outside the quantum dot, dielectric confinement and quantum confinement can be physically 
separated. In fact, by changing the dielectric environment far away from the dot, while 
keeping the same barrier material next to the dot, one can control and tailor the electronic 
properties (such as An,n+i and £|£p) without affecting quantum confinement (i.e. the single- 
particle energies and wave functions). 

We approximate the many-particle wave function of a. system of electrons in the 
conduction band of a quantum dot by a single Slater determinant constructed from the wave 
functions {ipi, i = 1 ■ ■ ■ N} of the N single-particle states occupied by the N electrons. The 
corresponding total energy is 
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Em = Eo + Y.{e, + Ef) n, + J^iAj - K,,,) n, n, , (8) 

i i<j 

where are the conduction-band single-particle energy levels, are the polarization self- 
energies, Jj j, Ki j are the electron-electron Coulomb and exchange energies, respectively, 
and Ui are the occupation numbers {J2i = N). The ground state \E'^ corresponds to the 
configuration that minimizes the total energy i?jv In Eq- ® we neglect: (i) the coupling 
between different Slater determinants (i.e. configuration- interaction effects), and (ii) the 
response of the single-particle wave functions i/ji to the electrostatic field (i.e. self-consistent 
effects). These assumptions are sufficiently accurate in small, three-dimensional structures 
in the strong-confinement regime [|16|-|TH]- 

The single-particle energies Si and wave functions ipi{r,a) are obtained here from the 
solution of the Schroedinger equation: 

[-V' + Vpsir)] ^i(r, a) = Si ^i{v, a). (9) 

The pseudopotential of the quantum dot V^^ij) is obtained from the superposition of screened 
atomic potentials, which are fitted [0 to reproduce the bulk experimental optical transition 
energies and effective masses, as well as the surface work function. Spin-orbit coupling is 
fully included in the solution of the Schroedinger equation. 
The interelectronic energies Jjj are given by: 

•^M = eE / |^.(r,a)|2$^.(r)dr, (10) 

where $j(r) is the electrostatic potential energy due to a charge distribution Pjij) = 
e X)<T |'^i(r, cr)P in a dielectrically inhomogeneous medium. $j(r) satisfies the Poisson equa- 
tion: 

V ■ e(r)V$j(r) = -Anepj{r), (11) 

where e(r) is the (position-dependent) macroscopic dielectric constant of the system. The 
Poisson equation is solved on a real-space grid using a finite-difference discretization of the 
gradient operator. The boundary conditions are obtained from a multipole expansion of the 
electrostatic potential fl^. The dielectric constant e(r) changes smoothly from ein to eout, 
with a transition region of the order of the interatomic bond-length. The interelectronic 
energy Jjj can be separated into two contributions: (a) the direct Coulomb energy Jf^j^^, 
which corresponds to the interaction between two electrons in the quantum dot as if the 
dielectric constant was uniform throughout the system, and identical to the macroscopic 
dielectric constant of the quantum dot; and (b) the polarization energy Jff which accounts 
for the effects of the dielectric discontinuity at the interface between the quantum dot and 
the surrounding material, and the ensuing surface polarization charge. 
The polarization self-energies are given by: 



-ipol 



^E./€(r,^)S(r)V'.(r,a)^ir, (12) 
E(r) = lim[G(r,r') - G'buik(r, r')], (13) 
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where G{r,r') is the Green's function associated with the Poisson equation [Eq. ([TT|)], and 
Gbuikli")!"') is the bulk Green's function. Here we use the analytical expression of S(r) for 
a spherical quantum dot [|1^] of dielectric constant embedded in a medium of dielectric 
constant eout- The singularity of S(r) at the surface of the dot is removed by applying a 
smoothing function 1 — e"'-'""^-' Z"" , where a is a broadening factor of the order of the bond 
length. 

We consider InAs spherical nanocrystals of diameter D = 30.3 and 42.2 A. The surface 



dangling bonds are passivated using a large-gap barrier material [|T^. Our analysis of the 
envelope functions extracted for the pseudopotential wave functions shows that the first 
electron level (el) is predominantly s-like, while the next 3 electron levels (e2, e3, and e4) 
are predominantly p-like. The first two hole levels {hi and h2) have an s-like envelope 
function, while the next two hole levels {h3 and M) have a p-like envelope function. Each 
single-particle energy level is doubly degenerate (because of time-reversal symmetry). 

The self-energies the polarization energies Jff, and the Coulomb energies Ji^"^^ 

of the 30.3 A diameter InAs nanocrystal are shown in Fig. 2 as a function of the external 
dielectric constant eout, for a few single-particle states i and j. We see that (i) both Ef°^ 
and Jff depend strongly on Cout, and vanish when eout = Cin (vertical arrows in Fig. 2). (ii) 
When eout > ^in the polarization energies Jff become negative, thus acting to diminish the 
electron-electron interaction, (iii) The dependence of and Jff on the identity of the 
orbitals i and j (e.g. s or p) is rather weak, as shown in the insets in Fig. 2. (iv) There is 
a critical value of eout (ecriticai ~ 4) such that for eout < Ccriticai the polarization energies Jff 
dominate over the Coulomb energies Jf^f^^- 

The charging energies fiN = Ej^ — -Eat-i, calculated from the total energies Ej^ given 
by Eq. (^, are shown in the central panel of Fig. 3 as a function of eout- The vertical 
arrow at the bottom of the figure denotes the value eout = fin, which divides the behavior 
into two domains: (i) In the weak screening regime (eout Cm) the charging energies are 
widely spaced, and their value depend strongly on eout- (ii) In the strong screening regime, 
on the other hand, the charging energies are closely spaced, and do not depend significantly 
on eout- The calculated charging spectrum is shown in Fig. 3 for eout = 1 (left-hand side) 
and eout = 20 (right-hand side), illustrating these two limiting behaviors. 

The electron and hole addition energies A7v,Ar+i (spacings between peaks in the charging 
spectra of Fig. 3) are given in Table I for a few values of eout- We see that (i) for a given value 
of eout, the addition energy for the third electron A2 3 is significantly larger than the addition 
energy for the second electron Af\. This can be explained by noting from Eqs. (^) and 

(I) that while aJ^J measures only the interelectronic repulsion, A2'^3 includes also the single- 
particle gap between the s-like state el and the p-like states e2, e3, and e4 {ee2—£ei = 0.40 eV 
for the D = 30.3 A nanocrystal and 0.36 eV for the D = 42.2 A nanocrystal). (ii) The 
addition energies of the remaining electrons (up to = 8) are approximately constant. The 
addition energy of the 9-th electron Ag^g is slightly larger, and reflects the single-particle 
gap between the jo-like shell and the next shell, (iii) While Coulomb blockade effects were 
often interpreted in terms of the "constant interaction" model, whereby the addition energies 
Aat^tv+i are constant and independent of A^, we see that in small nanocrystals deviations 
from this model are noticeable, and are due primarily to the existence of single-particle gaps 
comparable or larger than the interelectronic energies Jij. These deviations are magnified 
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when polarization effects are small (eout > Cin) • (iv) The electron and hole addition energies 
become smaller as the size D increases and/or the dielectric constant Cout decreases, (v) 
The quasi-particle gap e^^p depends strongly on eout, while the optical gap 6°^* does not. 
This is so because in the optical gap + Eg°' and J^fei tend to cancel [see Eq. (|^], so 
^gap ~ (^ei — — Jhi%- ('^i) The optical gap is smaller than the quasi-particle gap. For 
Cout ~ Cin the difference is Jhi^i- 

Figure 3 and Table I provide clear predictions of the charging energies and addition 
energies of InAs quantum dots. To compare with the experimental measurements of Banin 
et al. (Ref. [|10|)5 in which eout is an unknown quantity, we first fit our calculated A^°2 
the smaller dot with the experimental value Ai^2 = 0.22 eV, finding that eout = 6 gives a 
good fit (Table I). Using this value of eout, we then predict for D = 30.3 A (experimental 
data in parenthesis for D = 34A) e^Pp = 1.78 (1.75), A^"^^ = 0.24 (0.20), A^"^] = 0.23 
(0.22), aQ = 0.64 (0.71), and A^ = 0.24 (0.23). Using the same value of eout, our 
predictions for D = 42.2 A (experimental data in parenthesis for D = AAA) are: £|Pp = 1.38 
(1.38), aS"^^ = 0.16 (0.20), A^^ = 0.16 (0.17), Ag = 0.15 (0.14), A^^l = 0.51 (0.52), 

and A^l = 0.14 (0.14). We see that we have excellent agreement with experiment. Our 
theory can be further used to decompose the experimentally measured quantities into distinct 
physical contributions. For example, for D = 30.3 A the quasi-particle gap e^^^ = 1.78 eV 
includes [Eq. (|)] the single-particle gap 5ei — ^hi = 1-71 eV and the self-energy contribution 
"E^f + = 0.07 eV. The addition energy for the third electron A'^^^ = 0.64 eV includes 
[Eq. (H)] the single-particle contribution ee2 — Eei = 0.40 eV, the Coulomb contribution 
2^e°M2 - J?i,el = 0.17 eV, the polarization contribution 27^^% - JeiU = 0.07 eV, and a 
negligible self-energy contribution T,^2^ — Sei'- The exchange contribution i^ei,e2 is smaller 
than 0.02 eV, and can be neglected. 

In conclusion, we predict the effects of the dielectric environment on the electron and 
hole charging energies and on the addition spectrum of semiconductor quantum dots. We 
find that the charging energies and the addition energies depend sensitively on the dielectric 
constant eout of the surrounding material via the self-energies and the polarization 
energies Jff. When eout ^ Cin the charging energies are widely spaced in energy, and 
depend strongly on eout- When eout > Cin the charging energies are more closely spaced. Our 
calculations are in excellent agreement with recent spectroscopic results [0 for eout = 6. 

This work was supported by the U.S. DOE, OER-BES, Division of Materials Science, 
under Grant No. DE-AC36-98-GO10337. 
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TABLES 



TABLE I. Addition energies A7v,7v+i, quasi-particle gap £gip, and optical gap of InAs 

naiKHTvslals (iu cV) as a fuuciion of IIk^ diol(X'lric c(mslaul c,,^,. 



D = 30.3 A 


D = 42.2 A 


£gap = 1-71 eV 


£gap = 1-31 eV 



1 G 10 20 1 G 10 20 



Electrons 

A^*;^ 0.96 0.22 

A^^^ 1.45 0.64 

aJ^J 0.99 0.24 

aJ^J 0.98 0.23 

A^ 0.99 0.24 

A^^J 0.99 0.24 

Ag 0.99 0.24 

A^^^ 1.04 0.28 

A^o 1.00 0.25 

Holes 

A^^l 0.98 0.24 

A^"^] 0.97 0.23 

A^"^] 0.98 0.24 

A^^^ 1.02 0.23 

a'^^I 0.97 0.23 

Gaps 

e^Pp 2.37 1.78 

e°P* 1.56 



0.16 


0.11 


0.69 


0.57 


0.53 


1.05 


0.18 


0.13 


0.69 


0.17 


n 1 o 
U.lz 


0.7U 


0.18 


0.13 


0.69 


0.18 


0.13 


0.69 


0.18 


0.13 


0.69 


0.21 


0.17 




0.19 


0.14 




0.18 


0.13 


0.71 


0.17 


0.12 


0.71 


0.18 


0.13 


0.71 


0.16 


0.12 


0.73 


0.16 


0.12 


0.71 


1.72 


1.65 


1.84 


1.55 


1.54 





0.15 


0.10 


0.07 


0.51 


0.46 


0.43 


0.14 


0.10 


0.07 


U.io 


U.iU 


n HT 
U.Ui 


0.15 


0.10 


0.07 


0.14 


0.10 


0.06 


0.15 


0.10 


0.07 


0.16 


0.12 


0.09 


0.16 


0.12 


0.08 


0.16 


0.12 


0.09 


0.17 


0.13 


0.09 


0.16 


0.12 


0.08 


1.38 


1.32 


1.27 


1.22 


1.21 


1.20 
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FIGURES 




FIG. 1. Part (a) illustrates the process of loading three electrons into an otherwise neutral 
quantum dot. Part (b) shows the process of removing a single electron from a dot and placing it 
into another dot. 
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FIG. 2. The self-energy 11?° and the polarization energy J?° of an In As quantum dot (diameter 
D = 30.3 A) are shown as a function of the outside dielectric constant 6out- 

The insets show the 

differences SP°i - SP°i and Jf^^ - Jf°^ as a function of Cout- 

The vertical arrows indicate the value 

^out ~ ^in- 
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FIG. 3. The central panel shows the dependence of a the electron charging energies jii - ■ ■ /ig 
and of the hole charging energies //_!••• /I-q on the external dielectric constant eout- The vertical 
arrow indicates the value Cout = ^in- The side panels show the calculated charging spectrum in 
the case €out = 1 (left-hand panel) and eout = 20 (right-hand panel). The zero of the energy scale 
corresponds to the highest energy valence state. 
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